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<N ; ABSTRACT 

o 

We study the collimation of a highly magnetized jet by a surrounding cocoon 
that forms as a result of the interaction of the jet with the external medium. We 
show that in regions where the jet is well confined by the cocoon, current-driven 
^ ■ instabilities should develop over timescales shorter than the expansion time of 

the jet's head. We speculate that these instabilities would give rise to complete 
■ magnetic field destruction, whereby the jet undergoes a transition from high-to- 

low sigma above the collimation zone. Using this assumption, we construct a 
& ■ self-consistent model for the evolution of the jet-cocoon system in an ambient 

medium of arbitrary density profile. We apply the model to jet breakout in 
long GRBs, and show that the jet is highly collimated inside the envelope of the 

> ■ 

progenitor star, and is likely to remain confined well after breakout. We speculate 

C$ ■ 

that this strong confinement may provide a channel for magnetic field conversion 
in GRB outflows, whereby the hot, low-sigma jet section thereby produced is the 
source of the photospheric emission observed in many bursts. 

vb" 

(N 

in 

ON " 1. Introduction 

o; 

The relativistic outflows observed in many compact astrophysical systems are commonly 
thought to be powered by magnetic extraction of the rotational energy of a neutron star 
or an accreting black hole. The energy thereby extracted is transported outward in the 
form of Poynting flux, that on large enough scales is converted to kinetic energy flux. The 
mechanism by which magnetic energy is converted to kinetic energy has not been identified 
yet, but it is generally believed to involve gradual acceleration of the flow (e.g., Heyvaerts & 
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Norman 1989; Chiueh et al. 1991; Bogovalov 1995; Lyubarsky, 2009), impulsive acceleration 
(Granot et al. 2011; Lyutikov 2011; Granot 2012), and/or non-ideal MHD effects, specifically 
magnetic reconnection (Lyutikov & Blandford 2003; Giannios & Spruit 2007; Lyubarsky 
2010; McKinney & Uzdensky 2012). 

Astrophysical outflows appear to be highly collimated, and there have been many at- 
tempts to explain the observed collimation in different systems. Understanding the collima- 
tion process is important not only from an observational point of view, but also because in 
ideal MHD flows collimation and acceleration are intimately related (e.g., Begelman & Li 
1994; Vlahakis 2004). Magnetic fields can cause collimation via magnetic tension. However, 
collimation by magnetic tension alone is extremely slow (Eichler 1993; Begelman & Li 1994; 
Tomimatsu 1994; Beskin et al. 1998) and cannot account for the inferred collimation scales. 
Confinement by the pressure and inertia of an external medium has emerged as a promising 
alternative (Begelman 1995). The environments in which astrophysical jets propagate, e.g., 
accretion disk winds in the case of AGNs, or stellar envelopes in the case of long GRBs, are 
ideal for this purpose (e.g., Eichler 1983; Levinson & Eichler 2000; Bromberg & Levinson 
2007, 2009; Kohler et al. 2012). 

The effect of the external medium on the structure of MHD jets has been studied using 
semi-analytic models (Zakamska et al. 2008; Lyubarsky 2009, 2011; Kohler & Begelman 
2012) and numerical simulations (McKinney & Blandford 2009; Komissarov et al. 2007; 
Tchekhovskoy et al. 2010). However, these studies are restricted to steady state solutions in 
which the jet boundary is either treated as a rigid wall, or determined by external pressure 
with a prescribed profile. In more realistic situations it is expected that the jet will be sur- 
rounded by a hot cocoon that forms due to side flows of shocked matter from the jet's head, 
or a nose cone if magnetic pinching is important (Komissarov 1999). Indeed, numerical sim- 
ulations (Marti et al. 1997; Aloy et al. 1999; Hughes et al. 2002; Zhang et al. 2003; Morsony 
et al. 2007; Mizuta & Aloy 2009; Lazzati et al. 2009) and analytic models (Begelman & 
Cioffi 1989; Matzner 2003; Lazzati & Begelman 2005; Bromberg et al. 2011 (BNPS11)) of 
purely hydrodynamic jets indicate that under astrophysical conditions anticipated in GRBs, 
AGNs and microquasars, the surrounding cocoon significantly affects the structure and dy- 
namics of the jet. Attempts to simulate the propagation of a magnetized jet in an external 
medium have been limited to two-dimensional Newtonian jets (Clarke et al. 1986; Lind et al. 
1989), and relativistic jets with moderate magnetization, a ^ 1, under restricted conditions 
(Van Putten 1996; Komissarov 1999). Here a = B 2 /Anpc 2 , where B and p are the magnetic 
field strength and gas density in the jet. The results of such simulations should be treated 
with caution, as they cannot account for magnetic field dissipation in the shocked jet and in 
the nose cone that, as we will argue below, might be important. 
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In this paper we construct an analytic model for the propagation of a highly magnetized 
jet in an external medium. We suggest that a self-consistent treatment of the evolution of 
the jet-cocoon system may require a proper account of magnetic field dissipation at the jet's 
head. Such dissipation is implicitly invoked in our model. In section 2 we outline the basic 
model and its key features. In section 3 we consider the collimation of newly ejected jet 
material by the surrounding cocoon. In section 4 we compute the evolution of the jet-cocoon 
system under different conditions. In section 5 we consider the applications of our model to 
astrophysical systems. We conclude in section 6. 

2. The basic model 

Consider a magnetized jet propagating along the z-axis in a medium of density p a (z). 
We suppose that the jet is injected with a fixed opening angle Oo, a Lorentz factor jj, and a 
total power Lj. As the jet pushes its way forward, it accelerates the matter ahead of it and 
produces a strong forward shock. The jet is slowed down behind a reverse shock, to match 
the velocity of the shocked ambient matter at the contact discontinuity that separates the 
shocked jet plasma and the shocked ambient gas. In the case of a hydrodynamic jet the 
reverse shock is strong whenever the ambient density is large enough, such that the center of 
momentum frame moves at a Lorentz factor considerably smaller than jj. Then, the flow of 
matter into the jet's head through the forward and reverse shocks is balanced by a sideway 
flow that feeds a cocoon surrounding the jet. As long as the cocoon's pressure is sufficiently 
large, it deflects the streamlines of newly injected fluid and collimates the jet. As the system 
evolves, the cocoon expands and its pressure drops, until reaching a level at which it is too 
low to confine the jet. At which point in the course of the evolution this happens depends 
merely on the density profile of the ambient medium, the injected power Lj and the opening 
angle 8q at the injection point. 

The picture appears to be more involved in the case of magnetically dominated MHD 
jets. In the context of ideal MHD, the reverse shock is expected to be weak, implying that 
the velocity Uj = , jjf3j of the unshocked jet near the head should roughly match the velocity 
of the head, viz., Uj ~ Uh- Such a condition requires strong focusing of the section of the 
jet below the contact discontinuity, that can only be accomplished if the pressure of the 
confining medium roughly equals the pressure behind the forward shock. We shall propose 
in section I4TT1 that this might be accomplished through magnetic stresses, if some fraction of 
the magnetic field that enters the head through the reverse shock is advected into the cocoon 
and remains ordered over scales considerably larger than the cross sectional radius of the jet. 
An alternative possibility, discussed in section I4.2[ is that the section of the jet below the 



-4- 



contact discontinnuity is rendered unstable, thereby giving rise to rapid dissipation of the 
magnetic field there, and/or entrainment of shocked ambient matter, that can strengthen the 
reverse shock. The destruction of the magnetic field below the head may result from current 
driven-instabilities (Eichler 1993; Begelman 1998; Mizuno et al. 2009, 2012). Entrainment 
of matter may be driven by rapid growth of the Rayleigh- Taylor instability at the contact, 
which is expected when the head is accelerating (Lyubarsky 2010). The structure of the 
cocoon thereby formed consists of an outer part containing the shocked ambient plasma, 
and an inner part containing the lighter, shocked jet fluid (see figured]). 

As long as the pressure in the cocoon is sufficiently large it will collimate the jet. 
The transition from a freely expanding to a confined jet will occur at a radius at which 
the transverse ram pressure of the conical jet roughly equals the cocoon's pressure. If the 
transverse expansion of the unconfined jet is super-magnetosonic, then the streamlines of 
the jet will be deflected across a superfast tangential shock. Otherwise, the collimation will 
proceed smoothly through the formation of a rarefaction wave that propagates from the jet 
boundary inwards. The collimation of the jet by the cocoon is analyzed in section [3j The 
results derived there are used in section H] to determine the scales over which the cocoon 
significantly affects the evolution of the jet. 

The asymptotic structure of a relativistic, strongly magnetized jet has been analyzed 
recently by Lyubarsky (2009, 2011). He obtained analytic solutions of the transfield equation 
in the limit Qr ^> 1 (r is the cylindrical radius), that describe a rigidly rotating, steady jet 
confined by an external pressure having a power law profile p ex t{z) oc z~ K . He examined the 
behavior of the solutions in the regimes k < 2 and k > 2, and showed that when k < 2 the 
opening angle of the jet 9j decreases continuously such that the jet interior remains in causal 
contact (jjOj ~ 1) everywhere. As a consequence, the jet is accelerated and collimated until 
it roughly reaches equipartition, where a ~ 1. The jet's streamlines have a parabolic shape, 
r oc W 4 , with spatial oscillations superimposed on it if the jet is initially out of equilibrium. 
For k > 2 the jet becomes asymptotically radial, with the final opening angle 6joo depending 
solely on the pressure profile. In this case 7^00^00 > 1, so that the jet interior is not in causal 
contact. The asymptotic Lorentz factor is then practically limited to ;$ ^c/^joc/ 3 ' where 
cr denotes Michel's magnetization parameter. Since •yj 00 9joo > 1, the jet remains Poynting 
dominated in the far zone, viz., jj^ < a . As shown in Lyubarsky (2011), the above results 
hold not only for a cold jet, but also in the case of a magnetically dominated hot jet. 

The model outlined in the following sections assumes that above the transition zone, 
at radii where confinement by the cocoon has been communicated to the jet interior , the 
jet structure can be described by the solution derived in Lyubarsky (2009) for k < 2. This 
may be justified in the regime where the propagation of the jet's head is sub- or even mildly 
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relativistic, such that the vertical pressure distribution in the cocoon is sufficiently flat. As 
will be shown in section [3] below, the confined jet just upward of the transition region is 
out of equilibrium, and is expected to oscillate. These oscillations may eventually decay 
to an equilibrium state via production of shocks in the external gas flowing near the jet 
boundary. We suppose that near the position of the head the jet is in its equilibrium state, 
and use the equilibrium solution derived in Lyubarsky (2009) as a closure condition for the 
jet-cocoon system analyzed in section HI A schematic illustration of the model described 
above is presented in figure [TJ 



3. Transition to the collimation regime 

When the transverse ram pressure of the freely expanding jet roughly equals the cocoon's 
pressure, the streamlines of the jet will be deflected until the jet becomes confined by the 
gas in the cocoon. As mentioned above, whether the collimation of the jet proceed smoothly 
or through a formation of a tangential shock, depends on the transverse fast magnetosonic 
Mach number of the unconfined jet, Mf± = (v,j sin#o) / u f, where Uf is the fast magnetosonic 
4-speed, and 9q denotes the opening angle of the jet prior to its interaction with the cocoon. 
This is confirmed in appendix |Aj where the jump conditions of an oblique MHD shock are 
solved. 

For a magnetically dominated jet, the asymptotic Lorentz factor is limited to 7j £ (ctq / Oq) 1 ^ 3 
(Lyubarsky 2011), where do is the value of the magnetization parameter at injection (roughly 
the Michel's magnetization parameter). If the jet is injected cold, then jjaj = cr to 0{aj 1 ) 
(see Equation flA7|) ). yielding ■jjPjOo/ \fo~j ~ 1 in the asymptotic limit. Consequently, for 
a cold jet the transverse expansion is marginally sub- magnetosonic and we expect smooth 
collimation. The confinement of the jet in this case is communicated to the inner regions 
by a fast magnetosonic wave that propagates from the jet boundary inwards (Kohler & 
Begelman 2012). If the jet is injected hot then it may become superfast when interacting 
with the cocoon, and a tangential collimation shock will form. At the tangential contact 
discontinuity the pressure of the shocked jet layer must be equal to the cocoon's pressure p c . 
Then, assuming that the transverse momentum flux is roughly uniform across the shocked 
layer we obtain, using Equation (1A4|) . 



Pc = (hj + n;, J >] cos 2 ijij (1) 

for a superfast flow, viz., Mf± 3> 1- Here ipj is the angle between the fluid velocity and the 
shock normal (see, appendix lAl). 

Adopting p c = p (z/ Rl)~ k for the external pressure, Rl being the radius of the light 
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cylinder, and recalling that Lj = (hj + a j) pj^y 2 (3jC 3 7i tan 2 9qz 2 , yields 

cos tpj = Atan6 Q (z/R L y- K/2 , (2) 

with A = {iicpqR 2 l / Lj) 1 ! 2 ~ (Attpq/Bq) 1 / 2 , where Bq is the characteristic magnetic field at 
the light cylinder. We see that this condition is essentially the same as in the hydrodynamic 
case (Bromberg & Levinson 2009; Van Putten & Levinson 2012), as naively expected. In 
the small angle approximation cos ipj = r s /z — dr s /dz, where r s (z) is the cylindrical radius 
of the shock front. Substituting the latter expansion into Eq. fl2]) and solving the resultant 
ODE gives 

2A 

r s {z) = ztan6 - - z tan 6 [(z / R l ) 1 - k/2 - 1]. (3) 

2 — K 

The point z* at which the shock reaches the axis is determined from the condition r s (z = 
z*) = 0: 

2-K N V(l-«/2) 



" = M^r + 1 J ' (4) 

The latter result is used in section 14.21 to determine the scales over which collimation of the 
jet by the cocoon occurs (Equation 



If the motion of the head is slow enough, the cocoon is anticipated to be roughly isobaric. 
Then k — 0, and if A <C 1 we have z* / Rl ~ (Bq/Attpq) 1 ^ 2 . The radius of the jet at this 
distance is R* = Rl tan6 , o(-Bg/47rpo) 1 ^ 2 - The radius of an equilibrium cylinder is R e = 
Rl{BI/2-kpq) 1 /' 1 (Lyubarski 2009). Our analysis assumes that 7 J sin^ > \f a j/hj > 1- Since 
7j Rj/Rl-, it implies that sin^o > a/ aj/hj(R*/R L )^ 1 . At small angles sin^o ~ tan 6*0, 
and the above results yield R* £ ^J<jj/2hjR e > R e . If tan# < ^aj/hj(R e /Ri)~ 1 it means 
that the transverse expansion of the jet is sub-magnetosonic, so that collimation proceeds 
smoothly. 

Because Rj > R e at z*, we expect spatial oscillations of the jet to ensue above the 
transition region. To illustrate this, consider the propagation of a jet in a confining medium 
having a uniform pressure, p c = p = const. Suppose that at z = z the radius of the 
jet satisfies Rj = Ro and dRj/dz = 0. The solution of the transfield equation in this case 
reduces to (Lyubarsky 2009) 



Rj{z) = tfofcos 2 {^/3/2 A{z - z )} + {R e /R f sin 2 { ^j2A{z - ^ )}] 1/2 , (5) 

with A and R e as defined above. For R = R e the latter equation yields Rj(z) = R e = const. 
For any other values the jet oscillates. Such oscillations are present for any pressure profile 
with k < 2, when the initial jet radius deviates from the equilibrium value. These spatial 
oscillations may eventually decay and the jet radius will approach R e . 
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The above analysis applies as long as the Lorentz factor of the confined jet, jj ~ 
Re/Ri, is smaller than the maximum Lorentz factor achievable, jmax — 0"o- Once the 
cocoon's pressure drops to a level at which R e /R L > jmax, the jet will become kinetic 
dominated already during the collimation phase. The tangential shock, in that case, will 
become progressively stronger, converting bulk energy into kinetic pressure. In this regime, 
the confined jet will have low magnetization and a Lorentz factor jj ~ Oq 1 (BNPS11). 



4. Evolution of the jet-cocoon system 

The temporal evolution of the jet-cocoon system is determined by the density profile of 
the ambient medium and the parameters of the injected jet, assumed to be given. As ex- 
plained in section [2] above, matching of the jet and head Lorentz factors can be accomplished 
through magnetic pinching or, alternatively, non-ideal MHD effects that allow the formation 
of a strong reverse shock below the contact discontinuity. The evolution of the system in the 
former case is explored in section I4.1[ and in the latter case in section 14.21 

The energy momentum tensor of the unshocked jet can be expressed as 

Tf = (wj + fyuffl + ( Pj + tij/2)!r - ftjftf (6) 

where pj and Wj are the pressure and specific enthalpy, respectively, is the 4-velocity, 
and y/Airb^ = F* u u u is the magnetic field vector, F* being the dual electromagnetic tensor. 
Since well above the light cylinder the azimuthal magnetic field of the unshocked jet satisfies 
Bfj, ~ rQBp 3> Bp, where r is the cylindrical radius of the magnetic surface ^(r, z) and 
fi(^) its angular velocity, we can neglect the poloidal field B p . Then bj = (0, 0, bj, 0), where 
y/Anbj = Bj^/jj is the proper magnetic field. The jet power is obtained upon integration of 
the energy flux from the jet axis to its boundary, at r = Rj: 

Lj= / T° z 2nrdr= / ( Wj + tifirfcpfiirrdr. (7) 
Jo Jo 

We suppose that the structure of the confined jet well above the transition region can be 
approximated by the equilibrium solution obtained in Lyubarsky (2009) in the limit of rigid 
rotation, Q = const. Then, Lj ~ {Bq/ Att)ckR 2 l , where Bq = ^q/R\ is the characteristic 
magnetic field at the light cylinder, and 271-^0 = 2n^(Rj, z) is the total magnetic flux 
subtended by the jet. In terms of the cocoon's pressure at the jet boundary, p c (z), the 
cross-sectional radius of the confined jet is given by (Lyubarsky 2009) 
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and the Lorentz factor by jj ~ Rj/Ri- 

Let fih denotes the velocity of the contact discontinuity (the head) and 7^ the corre- 
sponding Lorentz factor. Momentum balance at the forward and reverse shocks gives 

/ \(w 3 + b 2 ) 1 2 1 2 {f3 J -(3 h ) 2 +p 3 + b 2 /2}27rrdr = / [w a7 2 /3 2 + p a }2irrdr } (9) 
Jo Jo 

where w a and p a are the enthalpy and pressure of the ambient medium, respectively, and Rf 

is the cross sectional radius of the forward shock. For simplicity, we shall consider a cold 

medium, p a = 0, and assume that the effective cross sectional radius of the forward shock 

roughly equals that of the jet. Then, Equation (jSJ) yields: — f3 h / ' f3j) 2 + p = 7^/^, in 

terms of the dimensionless parameters, 

/ = Lj/inRipae 3 ), (10) 
p=( Pj +b 2 /2)/p a c 2 . (11) 

Here, / represents the ratio between the total energy density of the jet, Lj/cnR 2 , and the 
rest-mass energy density of the surrounding matter, as measured in the rest frame of the 
ambient medium. Likewise, p is the ratio between the total jet pressure and the rest-mass 
energy density of the ambient matter. The solution for the head velocity reads: 

i-\i 2 -{i+p/mi-p-i)} 1 / 2 , 

I3 h = 1 i_p_\ — ■ ( 12 ) 

Note that for a highly magnetized jet, b 2 ^> Wj, one has, to a good approximation, p = 
l/{2^ 2 ). When the reverse shock is strong I ^> p and the latter solution simplifies to 

P h = 1 + /-V2' ( 13 ^ 

as derived in BNPS11. When the reverse shock is weak, specifically when r y 2 lf l {Pj— Ph) 2 < 1, 
Equation f|T2|) reduces to 

A - ff f (") 

In order to proceed, we need to specify the conditions in the inner cocoon. In what 
follows, we consider two different scenarios. 



4.1. Magnetized cocoon 



In this section we consider the possibility that the jet is pinched by magnetic stresses 
in the inner cocoon, assuming that a fraction of the magnetic flux that enters the head 
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through the reverse shock is advected into the inner cocoon. The magnetic flux entering the 
head is obtained upon contracting the tensor Q^ u = (b^Uj — bjUj) with the shock normal 
n u = j r (— (3 r , 0, 0, 1), measured in the frame of the ambient medium, where /3 r denotes the 
velocity of the reverse shock in the latter frame. This yields 

9f n v = B ilr % - r ) (15) 

in terms of the field Bj = ^jbj. The assumption that the return current is flowing in a thin 
sheet at some radius R c > Rj implies that the magnetic field in the inner cocoon scales as 
B c oc r _1 with cylindrical radius r. Thus, the rate at which magnetic flux in advected into 
the cocoon is roughly f3 r B c0 \ia(R c /Rj), where B c0 is the value of B c at the jet boundary, 
r = Rj. Flux conservation implies 

B c0 = £ s £ i7 r(&/& - VMRc/Rj)]' 1 - ( 16 ) 

In order not to crush the jet, the cocoon's pressure at the jet boundary, p c0 ~ B^ /8n, should 
not exceed the jet pressure, 5|/(87T7|). This yields the condition 

\n(R c /Rj) > £ B -Wr(P j/Pr - 1). (17) 



Suppose now that the separation of the forward and reverse shocks is much slower than 
the velocity of the head; that is, (3h — (3 r . Equation (IT?]) then implies £b(/3j — Ph) "C 1 
in the relativistic regime jj 3> 1. Unless the fraction £b is very small, the latter condition 
means that \(3j — (3h\ 1. To illustrate some key features of jet focusing by the inner 
cocoon, we compute the evolution of the system by invoking the extreme condition 7^ = jj. 
This, of course, is a gross approximation, as some velocity difference is required in order 
that magnetic flux will be advected into the inner cocoon, as indicated by Equation (I16p . 
Nonetheless, for (3j/(3h — 1 ~ {ijlh)^ 1 it may represent a reasonable approximation of a 
more realistic situation. Furthermore, we assume that the inner jet can be described by 
the equilibrium solution derived in Lyubarsky (2009), so that jj ~ Rj/R L . With (3j = (3h 
Equation (jUj) yields 7^ = p ~ l/^l], an d since = 7j we have, using Equation ffTOl . 

n = V* = {Jfa*T*" (18) 

for a density profile p a = p Zh~ a , where Zh = Zh/Ri- Comparing the latter with the un- 
magnetized case discussed in the next section, it is seen that the head velocity is larger by 
a factor 1.5zj/ 9 when magnetic pinching is effective. We remind the reader that the above 
results assume that the radius of the forward shock roughly equals the jet radius at the head. 
If magnetic pinching is effective mainly at a nose surrounding the head it would mean that 
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the jet is strongly focused only in the vicinity of the head. Then the assumption that the 
radius of the forward shock matches that of the head is questionable. 

Let us proceed by assuming that the outer cocoon is roughly isobaric, has a cylindrical 
geometry, and contains radiation-dominated gas, as described in detail in §4.21 Then, its 
pressure can be approximated by Equation (|22|) below. Using Equations ( ITS]) and (|22|) with 

Pc = Pc/ PaC 2 and Z?^ 1 — 1 ~ 1/27^, we arrive at 

Pc = (n/sy^p^li- 1 oc 4- (3+2Q)/3 . (19) 

Under the assumption B c oc r -1 , the pressure exerted on the jet by the inner cocoon is 
related to the thermal pressure p c of the outer cocoon, at r = R c , through 

Pco = Pc(Rc/Rj) 2 - (20) 

Now, for the equilibrium solution adopted here the jet radius Rj is given by Equation flH]) with 
p c replaced by p c0 , specifically Rj/R L = (2L j /7ccR 2 L p c0 ) 1/A = {2L j /TrcR 2 L p c ) 1/A (R c /R j ) l/2 . 
Combining the latter relation with Eqs. f[TB"j) and ( JT9"j) gives the radius of the current sheet 
in the inner cocoon: 

R c /Rj = 2(3/r?) 1 / 4 t 1 / 2 . (21) 

The above derivation implicitly assumes that in the vicinity of the head the relativistic 
jet is focused by magnetic pinching to a level at which the jet Lorentz factor can be con- 
tinuously matched to the head. One might naively conclude that if the ambient medium 
is sufficiently dense, then ultimately a bubble of submagnetosonic, high Poynting flux ma- 
terial will fill the expanding cocoon, quenching the inner jet to a narrow cylinder of radius 
Rj ~ Rl that propagates outwards at a subrelativistic speed, 7, ~ Rj/Ri — 1. In reality, 
such a structure is expected to be extremely unstable. What seems likely to happen is that 
current-driven instabilities will destroy the magnetic field in the inner jet before it is even 
deposited in the cocoon. Since the comoving growth time of the instability roughly equals 
the Alfven crossing time of the jet, t' ~ Rj/vA (Begelman 1998), the length scale over which 
the instability develops is Xqd ~ "fjcf ~ IjRj- If the jet is well collimated by the cocoon, 
then this scale is generally smaller than the distance Zh of the jet's head from the origin. It 
is, therefore, conceivable that the magnetization of the jet in some region below the head is 
strongly suppressed by the instability. We examine the consequences of such a process next. 

4.2. Unmagnetized cocoon 

If the magnetic field which is advected into the cocoon from the head dissipates on time 
scales much shorter than the expansion time of the head, then magnetic stresses in the inner 
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cocoon can be ignored. As shown below, in the absence of magnetic pinching the jet cannot 
be sufficiently focused, and its Lorentz factor may be much larger than the Lorentz factor 
7^ of the head, as in the pure hydrodynamic case. This means that the jet must decelerate 
across a strong reverse shock. We envision that the shock is formed in a low-sigma section of 
the jet, following the destruction of the incident magnetic field by the instabilities described 
above. 

Under the conditions envisaged the pressure in the cocoon is anticipated to be radiation- 
dominated, thus we adopt an adiabatic index of 4/3. Following BNPS11, we assume that the 
cocoon is roughly isobaric, and approximate its geometry as a cylinder of height Zh = c J (3hdt 
and cylindrical radius r c = cf (3 c dt, where f3 c = (p c / PaC 2 ) 1 ^ 2 is the lateral expansion velocity 
of the outer cocoon. The former assumption, that the energy distribution in the cocoon is 
approximately uniform, introduces a considerable simplification and may be justified when 
the motion of head is sufficiently slow. With the above approximations the cocoon's pressure 
is given by p c = E C /3V C , where E c = r/Lj f(l — f3h)dt is the total energy deposited in the 
cocoon and V c = 7rr 2 Zh its volume. The parameter 7] represents the fraction of the energy 
that enters the cocoon, as explained in BNPS11; at sufficiently low Lorentz factors of the 
head, 7^ < 2zh/Rj, for which it is in causal contact across its transverse direction, 77 = 1. 
Otherwise r\ = 2zh/( r yhRj)- Taking for simplicity r c = c J (3 c dt ~ (3 c ct, and likewise Zh = c(3ht, 
we obtain 

p = ^ (22) 
Pc 3ncR 2 L ffi ' 1 ] 

in terms of the fiducial time t = ctjR^. This expression is accurate up to an order unity 
factor that depends on the density profile of the ambient medium (see BNPS11). In the 
regime where the reverse shock is strong, (3h can be approximated by Eq. ffTBl . Substituting 
(3 C and j3h into the last equation and solving for p c yields 

1/2 

Pc 



Next, we suppose that near the head the cross-sectional radius of the jet can be approx- 
imated by Equation (JSj) and the Lorentz factor by jj = Rj/Ri- Solving Equations (IE]), ( jlQl) 
and (|23|) one finds 

2/3 



l = o. 2 Q V 2 / 9 1 -p— r 4 / 9 = l f h a /3 r 4 / 9 , 

\RlPa c J 



(24) 



where we adopt an ambient density profile of the form p a (z) = poz a , with zu = Zhj 'Rl, and 
define 
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The radius and Lorentz factor of the jet are found from (jSJ), (1231) and 

7 . = Rj /R L = (12/7 7 ) 1/6 Z 1/4 f L/3 = l.Sr/- 1 / 6 /; 74 ^ 6 ?/ 9 . (26) 

Result fl26|) holds at times at which 7 j < 7 max . In the regime where the latter condition 
is violated the jet becomes kinetic dominated already during the collimation process, as 
explained at the end of section [31 

The position of the head at time t is given by Zh = Pht- From Eq. (|T3|) it is readily seen 
that if / <C 1 then Zh ~ / 1//2 t, and if I 3> 1, Zh ~ t. In the former case we have 

I = / 9/7 4 6Q - 4)/7 , (27) 

and in the latter case (I 3> 1) 

i = /o4 6a " 4)/9 . (28) 

Consequently, for a > 2/3, Z increases with 2/, and the head accelerates. This is different than 
the hydrodynamic case where this happens at a > 2 (c.f., BNPS11). Note that when the 
head's motion becomes relativistic the Lorentz factor of the head scales as 7^ ~ / 1//4 / a/2 oc 
z^ a ~ 2 ^ 18 (see Eq. [13]). The pressure inside the cocoon satisfies p c cx z^ s with s = (2a + 8)/7 
when I <C 1 and s = (6a + 8)/9 when I ^> 1. From Eq. (jHJ) and the relation 7, = Rj/Rl 
we have 7 j cx 5 S//4 = ^( 3a + 4 )/ 18 for I 3> 1. The scaling of 7/1 and 7, with confirms that 
as long as the jet is confined by the cocoon it always accelerates faster than the head. 
The jet is considered collimated as long as z* < Zh, where z* is given by Eq. PJ with 
k = for an isobaric cocoon. The jet will be unconfined only when A < 1, for which 
z* ~ (Lj/tccRIpc) 1 ^ 2 cx z s J 2 . When s < 2 (a < 5/3) the head advances faster than z*, and 
the jet becomes confined at distances z > z c , where 

z c = (l.Q v - 2 /% /2 ) 9 ^' 3a \ (29) 

When s > 2 (a > 5/3), the jet is confined at z < z c , and becomes unconfined at z > z c . 

The model presented above implicitly assumes that current-driven instabilities lead to 
magnetic field dissipation above the transition zone, at z > z*. This assumption is justi- 
fied provided the instability growth time is shorter than the expansion time of the head. 
As argued at the end of section 14.11 the instability growth length is Xcd ~ IjRj^A — 
R 2 /(PaRl), where [3a is the Alfven speed in units of c, for which Equation (JS]) yields 
Xcd/Rl — (ZLj/ircRlpc) 1 / 2 ^ 1 = s/^Pj^z*. Now, as long as the jet is confined by the 
cocoon z* < Zh, hence the magnetic field in the jet has sufficient time to dissipate before 
it reaches the head provided Pa > z*/zh- In particular, in the regime where the head is 
sub-relativistic we find Xcd/zh ~ P^t^ 3 ~ Pj^^^^h 3 ^ ? ' anc ^ m the re gi me where the 
head is transrelativistic Xcd/zh ~ P^Iq • 



13 



It is worth noting that in a relativistically hot, pure hydrodynamic flow the cross- 
sectional radius and Lorentz factor scale as jj oc Rj oc pj 1 ^ with pressure pj. By comparing 
with Equation it is seen that conversion of magnetic energy to kinetic energy in the 
confinement region does not change the scaling of the outflow parameters, so that the use 
of Equation ([8]) is justified even if the jet becomes kinetic dominated in the vicinity of the 
head. 



5. Applications 

We examine first the application of the above results to GRBs, assuming an unmag- 
netized cocoon. In the context of the collapsar scenario for long GRBs, the jet propagates 
inside the envelope of a massive star before breaking out to produce the observed signal. For 
illustration, we invoke a WR star of mass M ~ 1OM and radius R+ ~ Rq. The density 
profile in the stellar envelope may be expressed as p a (z) = p(z/R*)~ a 1 with a < 3, where 
the average density is roughly p ~ 5(3 — a) (M / 10 M Q )(R*/R Q )~ 3 g cm -3 . From Equation 
fl25|) we obtain 



/o = 2 x ioy/ 9 



2a/3 



R 2 

L7 



2/3 



(3 - a)M 
10 M 



-2/3 



Uo 



(30) 



where R L = 10 7 R L7 cm. For a = 2 and the above choice of M, R± and Rl we have, using 
(EZD with Tj = l, 



I ~ 5 x 10" 



L J52 
R L7 



6/7 



=8/7 
z h 



123 



52 



^L7 



6/7 



8/7 



(31) 



in the sub-relativistic regime, and from (T28]) 



Z ~ 1.5 x 10- 2 r] 2/g 



52 



^L7 



2/3 



-8/9 



40r/ 2 / 9 



f^52 
^L7 



2/3 



Zh 

R* 



8/9 



(32) 



in the relativistic regime (I ^> 1). The Lorentz factor of the head can be expressed as 7^ ~ 
2(Lj 52 / R^) 1 ^ 6 ( z h/ R*) 2 ^ 9 ; and it is seen that for values of the jet power inferred from observa- 
tions, Lj52 ~ 1, the motion of the head is sub-to-mildly relativistic inside the stellar envelope. 
For the above choice of parameters the head becomes relativistic at Zh > Zt — 0.0 14R*. The 
jet Lorentz factor at the location of the head is given by jj ~ 40(L J - 5 2/-R| 7 ) 1//14 (z/ l /-R^) 3 / 7 at 
z h < zt, and jj ~ Q2(L j52 /R 2 L7 ) 1/6 (z h /R^) 5/9 at z h > z t . From Eq. ([29]) we also have 



Zc/R* 



70 



L j52 
R L7 



-3 



(33) 



- 14 - 

Thus, the jet will be collimated all the way to the edge of the stellar envelope provided 
Lj 52 /R 2 L7 < 4. In fact, jets of sufficiently low power may remain confined well above the edge 
of the stellar envelope by the surrounding matter that breaks out of the star with the jet. 
This raises the question of whether the deconfinement of the jet at breakout is sudden enough 
to lead to re-acceleration to Tj/Oj ^> 1, as proposed recently (Tchekhovskoy et al. 2010; 
Komissarov et al. 2010). The dissipation of the magnetic field in the confinement region 
may persist after breakout, until the confinement relaxes to the point that current-driven 
instabilities no longer have time to operate. From the above results we obtain \cdI z h ~ 
0.55 / 9^ 1 (Lj52/-R| 7 ) 1 / 3 (2;/ l /-R*) 1//9 in the transrelativistic regime (7^ > 1). Thus, under the 
assumption that the cocoon is isobaric it seems that the instability is marginal near and 
beyond the edge of the star, and it is unclear whether it will indeed become destructive 
there. On the other hand, a non-uniform pressure distribution in the cocoon would lead 
to additional focusing of the jet above the collimation zone and, hence, a shorter growth 
length of the instability. The hot, low-sigma outflow produced by the destruction of the 
magnetic field well inside the star and, conceivably, also after breakout will eventually reach 
the photosphere, and may be the source of the photospheric emission observed during the 
prompt phase. 

Blazar jets propagating in a medium of constant density, n = Po/ m P (measured in cgs 
units), will be collimated at distances z > z c , with 

z c ^ 10 ( L p 2 5 ) pc, 

where we adopt R L = 10 14 R L i4 cm. For n ~ 1 cur 3 , z c is much larger than the scale over 
which the magnetic field dissipates, as inferred from observations. The cocoons observed on 
tens of kpc scales form, most likely, in the pure hydrodynamic regime (BNPS11). However, 
the average density in the vicinity of the broad line region may be larger by several orders 
of magnitudes than the average ISM density, so that collimation by the cocoon may occur 
on much smaller scales, shortly after the activation of the central engine. We note that 
the Lorentz factor of the confined jet at z > z c is 7-,- ~ 67(Lj^/ Rj jU no) 3 ^ w (zh/ z c ) 2 ^ 9 . The 
latter exceeds the characteristic values inferred from observations, 7^ ~ 10 — 50, implying 
that conversion of magnetic energy to kinetic energy can occur by collimation alone, even on 
sub-parsec scales if the ambient density in the vicinity of the central engine is large enough, 
viz., n > 10 3 . 



(34) 
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6. Conclusion 

We have constructed an analytic model for the propagation of a magnetically dominated 
jet in an external medium having an arbitrary density profile, and examined the conditions 
under which the jet can be collimated by the cocoon surrounding it. Our model assumes that 
in the region where the jet is well-confined by the cocoon, its structure can be described by 
the equilibrium solution derived in Lyubarsky (2009) in the equilibrium regime. This solution 
is employed as a closure condition for the jet-cocoon system of equations. We analyzed two 
different evolutionary tracks: 

The first one, outlined in §4.1[ assumes that some fraction of the magnetic field that 
enters the jet's head through the reverse shock is advected into the cocoon and remains 
ordered. The inner jet is focused by magnetic pinching to a level at which the jet Lorentz 
factor can be matched to the head through a weak reverse shock. Under these assumptions 
a self-consistent solution of the jet-cocoon equations can be obtained. This solution may 
correspond to the nose cone revealed in two-dimensional numerical simulations of moderately 
magnetized jets (Komissarov 1999). However, this scenario requires that the system maintain 
a high degree of axisymmetry over times longer than the expansion time of the head, which 
is questionable. We argued, at the end of §4.1[ that for a dense enough medium such a 
configuration should be unstable, and should lead to a rapid dissipation of the magnetic 
field in the inner jet and the surrounding cocoon. 

In the second one, which we find more likely, the magnetic flux in the confined jet 
dissipates via current- driven instabilities before it reaches the head, and the jet undergoes a 
transition from high-to-low sigma above the collimation zone. A strong reverse shock forms 
in the low-sigma section of the jet, allowing matching of the jet and head Lorentz factors, 
as in the pure hydrodynamic case. Since the side flow that feeds the cocoon is weakly 
magnetized, magnetic pinching is unimportant. Instead, the jet is confined by the pressure 
of the gas contained in the inner cocoon. Collimation commences at a radius at which the 
transverse ram pressure of the unconfined jet roughly equals the cocoon's pressure, followed 
by spatial oscillations of the confined jet that decay to the equilibrium state by dissipative 
processes. The collimation proceeds smoothly if the transverse expansion of the unconfined 
jet is sub-magnetosonic, or through formation of a (weak) superfast tangential shock if the 
transverse expansion is super-magnetosonic. At the onset of collimation the jet may still be 
highly magnetized (a ^> 1), unlike the pure hydrodynamic case considered in BNPS11. This 
leads to scalings different than those derived by BNPS11. For example, for transrelativistic 
propagation in an ambient density profile p a oc z~ a , we find that the head Lorentz factor 
evolves according to 7^ oc z^ a 2 ^ 18 in the high-sigma case, versus 7^ oc 2 ^ 10 in the pure 
hydrodynamic case. 



- 16 - 



As long as the jet is well-confined by the cocoon, the growth time of the current-driven 
instabilities is shorter than the expansion time of the head. As a consequence, complete 
destruction of the magnetic field is expected in the confinement region below the head. The 
hot, low-sigma matter thereby produced may be a source of high-energy radiation when 
approaching the photosphere. In the collapsar scenario for long GRBs, we find that the jet 
will remain well-confined throughout its propagation in the envelope of the progenitor star, 
and perhaps even well above the envelope. For a reasonable density profile, the criterion for 
the growth of the instability is found to be marginal near the edge of the envelope, so further 
analysis is required to quantify the likelihood that the instability will become disruptive. If 
it does, then magnetic field dissipation may persist for times longer than the duration of the 
breakout phase. At any rate, the hot, low-sigma matter produced inside and, conceivably, 
above the stellar envelope will eventually reach the photosphere and radiate. The photo- 
spheric emission observed in the prompt phase of many bursts may be a signature of this 
mechanism. Since magnetic field dissipation commences well inside the envelope, at large 
optical depths, there is sufficient time to generate the radiation entropy required to explain 
the sub-MeV peaks; typically, a Thomson depth of r > 10 3 is required for complete thermal- 
ization (Levinson 2012). The overall shape of the spectrum emitted from the photosphere 
would depend on the dissipation profile below the photosphere (Levinson 2012; Beloborodov 
2012), which in turn depends on the density profile of the progenitor star and other details. 
Further dissipation of the bulk energy of the weakly magnetized fluid near the photosphere 
may occur via formation of internal or collimation shocks. 

AL acknowledges support from an ISF grant for the Israeli center for high energy astro- 
physics, and thanks the Fellows of JILA for their hospitality during a sabbatical visit. MCB 
acknowledges support from NSF grant AST-0907872 and NASA Astrophysics Theory grant 
NNX09AG02G. 



A. Oblique MHD shocks in super-fast flows 

The energy momentum tensor of the upstream flow, as measured in the shock frame, is 
expressed as 

TT = (w 1 + 6?K< + (px + 6?/2)(T - « (Al) 

where is defined below equation §6§. We consider a planar shock and choose our coordinate 
system such that the velocity of the unshocked flow is given by (3 = (f3 x , 0, (3 Z ), and the shock 
normal by n = (n x , 0, n z ). For simplicity we assume that magnetic field of the unshocked flow 
just upstream of the shock is purely toroidal. Then b± = (0, 0, &i, 0), where V^brfei = B^/'ji 
is the proper magnetic field. In terms of the angle ip between the jet velocity and the shock 
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normal (cos^ = n • /3), the jump conditions are written as 



Pi7i/3i cos^i = P2I2P2 cos ip 2 , (A2) 

(hi + (Ti)pi7i/3i COS ^1 = (h 2 + 02)P2llfi2 cos ^ 2 (A3) 
(/l! + CTl)pi7i^ COS 2 -01 + Pi + Plfl/2 = (/l 2 + <T 2 )p272/ 3 2 COs2 ^2 + P2 + P20"2/2 (A4) 

/?i sin *0! = /3 2 sin ip 2 (A5) 
0"i/pi = cr 2 /p 2 . (A6) 

Here subscript 2 refers to shocked fluid quantities, a = 6 2 /p is the magnetization, and 
h = w/p is the enthalpy per baryon. Note that y/ojh is the Alfven 4- velocity. Equations 
( 1A2j) and (IA3I) can be combined to yield 

(/ii + (Ti)7i = (/i 2 + ^2)72- (A7) 

After some algebraic manipulations, the jump conditions can be reduced to a cubic equation 
for the variable x = u^ ± = 7|/3| cos 2 tp2- 



a^x + a 2 x + a\x + ao = 0, (A8) 



with 



a 3 = 16[c 2 -(l-/3 1 2 sin 2 ^ 1 )], (A9) 

a 2 = 16c 2 - 8cic 2 /3icos^i - 8(1 - /3 2 sin 2 ^i), (A10) 

ai = (? 2 $\ cos 2 ipi — 8cic 2 /3i cos-^i — (1 — 01 sin 2 ipi), (All) 

a = C 2/3 2 cos 2 ^i, (A12) 

where C\ = fii cos^x[l + (pi/(h\ + a\) + c 2 /2)/(w 2 cos 2 ^x)] and c 2 = cf\/(gi + Solutions 
of (1A8[) obtained numerically for /i x = 4p 1 /p 1 , cos-^i = 0.1, = 10, and different 

values of the fast magnetosonic Mach number, Mf± = Ui±/uf, Uf — a/ {h\ + 3a\)/2hi § , 
are exhibited in fig El As expected, shock solutions exist only for Mf± > 1. In the limit 
7i >> 1, M 2 _|_ >> 1 Eq. flA8[) can be solved analytically to yield 

2 8a 2 + lOfeig + fr 2 [64a 2 + a) 2 + 20(x(fr 1 + a) + Z^] 1 / 2 

M2± " le/i^a + ZiO le/n^ + Zix) ' 1 j 



1 In general, the fast magnetosonic speed c/ satisfies = (jp + b 2 )/(ph + b 2 ), where 7 is the adiabatic 
1^ = c^/(l - cj) = {h + 3a)/2h. 



index. In the limit h — 4p/p, 7 = 4/3, it reduces to A = (h/3 + cr)/{h + a), from which we obtain 
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which in the special case cos^i = 1, hi — 1 reduces to that obtained by Kennel and 
Coroniti (1984), and in the limit <j\ — 0, h\ — 1 reduces to U2± = 1/ \/8- Note that 
72 = (1 + v>2±)/0- ~ Pi sin 2 ^l)- As pointed out by Komissarov (2012), the proper condition 
for the approximation flA13[) is Mf± » 1, not just u±± » 1. In the limit <j\ » hi the 
solution (1A13j) simplifies to 



«a± = Mhi), (A14) 

; \ 1/2 

o\ + h\ sin 2 ifji 

P2 lifii , 
cos ip i, 



sin -02 = ( — — — — 2 ; ) sinipx, (A15) 



Pi \J(J\lh x 



Ptot =P2 + b 2 2 /2 = \U ^f 1 + 1 ] Mi7iA 2 cos^ 2 , (A17) 

in agreement with the high Mach number limit of figure [2j In fact, figure [2] indicates that 
the latter solution is a good approximation already at modest Mach numbers, Mf± = a few. 
The deflection angle of streamlines across the shock, 8 — ip2 — ipi, is readily obtained from 
the above: 

sinS = ( Vl + friM- M sin2 ^ ~ I(/ ll / (Tl ) 3 / 2 sin2^ 1 . (A18) 
\2y/<T 1 /h l + sin 2 </>i / 4 

The shock compression ration is r = Pi/ P2 — Piy/1 + h\/o\ y and it is seen that the in the 
limit <J\/h\ » 1 the shock is always weak. We emphasize that Eqs.( lAl~4l) - (1A17 l ) hold not 
only in the case of hot upstream flow but for any h\. 
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Fig. 1. — Schematic illustration of the jet-cocoon system. 



-22- 



1.2 









■ i 




i 




\ 




- \ 




\ 




\ 




\ 




\ 




\ 




\ 




\ 












; s v v pi >P2 












_ • — ^ 

" ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ 


, , I < < < I < < 



2 4 6 8 10 

Fig. 2. — Normal component of the downstream 4-velocity normalized to ua = \fo\fh\ 
(solid line), and proper density ratio p\j p 2 (dashed line), versus fast magnetosonic Mach 
number Mf±, for a hot upstream flow (hi = Api) with ai/hi = 10 and incidence angle 
(measured with respect to the shock normal) cos-^i = 0.1. 



